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Abstract. The paper is concerned with the boundary controllability of entropy weak solutions to 
hyperbolic systems of conservation laws. We prove a general result on the asymptotic stabilization 
of a system near a constant state. On the other hand, we give an example showing that exact 
controllability in finite time cannot be achieved, in general. 







1 - Introduction 



Consider an n x n system of conservation laws on a bounded interval: 

u t + f{u) x = t>0, xe]a,b[. (1.1) 

The system is assumed to be strictly hyperbolic, each characteristic field being either linearly 
degenerate or genuinely nonlinear in the sense of Lax [8] . We shall also assume that all characteristic 
speeds are bounded away from zero. More precisely, let / : ft i— » M n be a smooth map, defined on 
an open set f£ C M n . For each u € f2, call X\(u) < ■ ■ ■ < X n (u) the eigenvalues of the Jacobian 
matrix Df(u). We assume that there exists a minimum speed c > and an integer p € {1, . . . , n) 
such that 

\i(u) < if i < p, , 12 , 

\i(u) > if % > p, ^ ' ' 

\\{u)\ > c > uefl. (1.3) 

By (1.2), for a solution defined on the strip t > 0, x G]ct, b[, there will be n — p characteristics 
entering at the boundary point x = a, and p characteristics entering at x = b. The initial- 
boundary value problem is thus well posed if we prescribe n — p scalar conditions at x = a and p 
scalar conditions at x = b [11]. See also [1, 2] for the case of general entropy- weak solutions taking 
values in the space BV of functions with bounded variation. 

In the present paper we study the effect of boundary conditions on the solution of (1.1) from 
the point of view of control theory. Namely, given an initial condition 

u(0,x) = (f)(x) x e]a,b[ (1.4) 

with small total variation, we regard the boundary data as control functions, and study the family 
of configurations 

K(T) = |n(T,-)} Cl/Qa^]; M n ) (1.5) 
which can be reached by the system at a given time T > 0. 

Beginning with the simplest case, consider a strictly hyperbolic system with constant coeffi- 
cients: 

u t + Au x = 0, (1.6) 
where A is a n x n constant matrix, with real distinct eigenvalues 

Ai < • • • < A p < < A p+ i < ■ • • < A n . 

Call 

b — a 
t = max — — p- 

* |Aj| 

the maximum time taken by waves to cross the interval [a, b]. In this case, it is easy to see that 
the reachable set in (1.5) is the entire space: 1Z(T) = L 1 for all T > r. In other words, the 
system is completely controllable after time r. Indeed, for any T > r and initial and terminal data 
(f), tp G L 1 ([a, b]; M n ) , one can always find a solution of (1.4), defined on the rectangle [0, T] x [a, b] 
such that 

u(0,x) = <p(x), u(T,x) = ip(x) x€[a,b]. 
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Such solution can be constructed as follows. Let l\, . . . , l n and r\, . . . , r n be dual bases of right and 
left eigenvectors of A so that k ■ rj = 5ij. For i = 1, . . . , n, let Ui(t, x) be a solution to the scalar 
Cauchy problem 

u i)t + XiUi !X = 0, 

{k ■ 4>{x) if x G [a, 6], 

h-iP(x + \ t T) if s G [a -XiT, b- XiT], 

otherwise. 



Then the restriction of 



u(t,x) = ^2uj(t,x)rj 



to the interval [0,T] x [a, 5] satisfies (1.6) and takes the required initial and terminal values. Of 
course, this corresponds to the solution of an initial-boundary value problem, determined by the 
n boundary conditions 

{U ■ u(t, a) = Ui(t, a) i=p+l,...,n, 
k -u(t,b) =Ui(t,b) i = l,...,p. 

This result on exact boundary controllability has been extended in [9, 10] to the case of general 
quasilinear systems of the form 

u t + A(u)u x = 0. 

In this case, the existence of a solution taking the prescribed initial and terminal values is obtained 
for all sufficiently small data (p, ip € C 1 . 

Aim of the present paper is to study analogous controllability properties within the context of 
entropy weak solutions 1 1— ► u(t, •) € BV. For the definitions and basic properties of weak solutions 
we refer to [4]. For general nonlinear systems, it is clear that a complete controllability result 
within the space BV cannot hold. Indeed, already for a scalar conservation law, it was proved in 
[3] that the profiles tp € BV which can be attained at a fixed time T > are only those which 
satisfy the Oleinik-type conditions 

nx) - (x-a)f»(4>(x)) f ° ra - e - * eM - 

For general n x n systems, a complete characterization of the reachable set 1Z{T) does not 
seem possible, due to the complexity of repeated wave-front interactions. 

Our first result is concerned with stabilization near a constant state. Assuming that all 
characteristic speeds are bounded away from zero, we show that the system can be asymptotically 
stabilized to any state u* G with quadratic rate of convergence. 

Theorem 1. Let K be a compact, connected subset of the open domain Q C M n . Then there exist 
constants Cq,5,k > such that the following holds. For every constant state u* G K and every 
initial data u(0) = <j> : [a, b] i-> K with Tot.Var.{4>} < 5, there exists an entropy weak solution 
u = u(t,x) of (1.1) such that, for all t > 0, 

Tot.Var.{u(t)} < C e" 2 ^ , (1.7) 
\\u(t,x)-u*\\ TOO <C e- 2Kt . (1.8) 
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The proof will be given in Section 2. An interesting question is whether the constant state 
u* can be exactly reached, in a finite time T. By the results in [9], this is indeed the case if the 
initial data has small C 1 norm. On the contrary, in the final part of this paper, we show that exact 
controllability in finite time cannot be attained in general, if the initial data is only assumed to be 
small in BV . 

Our counterexample is concerned with a class of strictly hyperbolic, genuinely nonlinear 2x2 
systems of the form (1.1). More precisely, we assume 

(H) The eigenvalues \i(u) of the Jacobian matrix A(u) = Df(u) satisfy 

-A* < Ai(tt) < -A* < < A* < A 2 («) < A* . (1.9) 
Moreover, the right eigenvectors ri(u), r 2 {u) satisfy the inequalities 

£>Ai-n>0, L>A 2 -r 2 >0, (1.10) 

nAr 2 <0, n A (Dn • n) < 0, r 2 A (Dr 2 ■ r 2 ) < 0. (1.11) 

A partucular system which satisfies the above assumptions is the one studied by DiPerna [7] : 

Pt + {up) x = 0, 
u 2 x 



u t+[^ + TP 1 ' 1 ) =0, 
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with 1 < 7 < 3. Here p > and u denote the density and the velocity of a gas, respectively. 

The last two inequalities in (1.11) imply that the rarefaction curves (i.e. the integral curves of 
the vector fields r\,r 2 ) in the (u\,u 2 ) plane turn clockwise (fig. 1). In such case, the interaction of 
two shocks of the same family generates a shock in the other family. 




Theorem 2. Consider a 2 x 2 system satisfying the assumption (H). Then there exist initial data 
4> : [a, b] I— > iR 2 having arbitrarily small total bounded variation for which the following holds. For 
every entropy weak solution u of (1.1), (1-4), with Tot.Var.{u(t, •)} remaining small for all t, the 
set of shocks in u(t, •) is dense on [a, b], for each t > 0. In particular, u(t, •) cannot be constant. 

As a preliminary, in Section 3 we establish an Oleinik-type estimate on the decay of positive 
waves. This bound is of independent interest, and sharpens the results in [5], for systems satisfying 
the additional conditions (H). 

As a consequence, this implies that positive waves are "weak", and cannot completely cancel 
a shock within finite time. The proof of Theorem 2 is then achieved by an induction argument. We 
show that, if the set of 1-shocks is dense on [0, T] x [a, b], then the set of points Pj = (tj,Xj) where 
two 1-shocks interact and create a new 2-shock is also dense on the same domain. Therefore, new 
shocks are constantly generated, and the solution can never be reduced to a constant. Details of 
the proof will be given in Section 4. 

As in [9], all of the above results refer to the case where total control on the boundary values 
is available. As a consequence, the problem is reduced to proving the existence (or nonexistence) 
of an entropy weak solution defined on the open strip t > 0, x G]o, b[, satisfying the required 
conditions. This is a first step toward the analysis of more general controllability problems, where 
the control acts only on some of the boundary conditions. We thus leave open the case where a 
subset of indices I C {1, . . . , n} is given, and one requires 

li -u(t,a) = f { i = p+l,...,n, 

l-. u (t b) = i a ^ t) if i = l » 

for some control functions on acting only on the components i G I. 



Throughout the following, we denote by ^(u), U{u) the right and left i-eigenvectors of the 
Jacobian matrix A(u) = Df(u). As in [4], we write a i— ► Ri(a)(u ) for the parametrized i- 
rarefaction curve through the state Uq, so that 

^-R i (o)=r i (R i {o-)), Ri(0) = u o . 

The i-shock curve through uq is denoted by a i— ► Si(a)(uo). It satisfies the Rankine-Hugoniot 
equations 

f(Si(a)) - f(u ) = Ai(cr) (Si(a)-uo) 

for some shock speed A^. We recall (see [4], Chapter 5) that the general Riemann problem is solved 
in terms of the composite curves 

%(uo)(a) = { *i U °)W> ^ a -°' (1-12) 
n UA 1 \ Si(u )(a), if a < 0. v ' 
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2 - Proof of Theorem 1 



The proof relies on the two following two lemmas. 

Lemma 1. In the setting of Theorem 1, there exists a time T > such that the following holds. 
For every pair of states u>,u' £ K there exists an entropic solution u = u(t,x) of (1.1) such that 

u(0,x)=u, u(T,x)=lo' for all x € [a, b}. (2.1) 

Proof. Consider the function 

$(<7i, . . . ,a n ; v,v') = ^ n (o- n ) o ■ • • o V p+1 (a p+1 )(v') - ^ p (a p ) o • ■ ■ o *i(<7i)(t;). (2.2) 

Observe that, whenever v = v' , the n x n Jacobian matrix dQ/dai ■ ■ ■ a n computed at o\ = 
o~2 = ■ ■ ■ = o- n = has full rank. Indeed, the columns of this matrix are given by the linearly 
independent vectors —ri(v),...,—r p (v), r p +i(f), . . . , r n (v). By the Implicit Function Theorem 
and a compactness argument we can find 6 > such that the following holds. For every v, v' € K, 
with \v — v'\ < 5, there exist unique values <ti, . . . , a n such that 

v" = ^ n {a n ) o • • • o m p+l {a p+l ){v') = $ P K) o • • • o *i(<7i)(t;) . (2.3) 

Defining the time 

b — a . . 

t = max sup . . , (2.4) 

l<i<n uen \Ai{U)\ 

we claim that there exists an entropy weak solution u : [0, 2r] x [a, b] i— > such that 

u(0,x) = v, u(2r,x) = w / . (2.5) 




figure 2 
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The function u is constructed as follows (fig. 2). For t € [0,r] we let u be the solution of the 
Riemann problem 

»<°'*H:» I 111 < 26 > 

Moreover, for t € [r, 2r], we define u as the solution of the Riemann problem 

It is now clear that the restriction of u to the domain [0, 2r] x [a, 6] satisfies the conditions (2.5). 
Indeed, by (2.3), on [0, r] the solution u contains only waves of families < p, originating at the 
point (0,b). By (2.4) these waves cross the whole interval [a, b] and exit from the boundary point 
a before time r. Hence u(t,x) = v" . Similarly, still by (2.3), for t G [r, 2r] the function u contains 
only waves of families > p+ 1, originating at the point (r, a). By (2.4) these waves cross the whole 
interval [a, b] and exit from the boundary point b before time 2r. Hence u(2t, x) = v'. 

Next, given any two states u,ui' £ K, by the connectedness assumption we can find a chain of 
points loq = uj,u>i, . . . ,Un = oj' in K such that \oJi~ a>j_i| < 5 for every i = 1, . . . , N. Repeating the 
previous construction in connection with each pair of states (uJi-i, tOi), we thus obtain an entropy 
weak solution u : [0, 2Nt] x [a, b] i— > Q that satisfies the conclusion of the lemma, with T = 2Nt. 

□ 



In the following, we shall construct the desired solution u = u(t, x) as limit of a sequence of 
front tracking approximations. Roughly speaking, an e-approximate front tracking solution is a 
piecewise constant function n £ , having jumps along a finite set of straight lines in the t-x plane say 
x = x a (t), which approximately satisfies the Rankine-Hugoniot equations: 



^2 \ f(u(t,x a +)) - f(u(t,x a -)) - x a (u(t,x a +) - u(t,x a -)) 

a 

for all t > 0. For details, see [4], p. 125. 



< e 



Lemma 2. In the setting of Theorem 1, for every state u* E Vt there exist constants C,8q > 
for which the following holds. For any e > and every piecewise constant function u : [a, b] i— ► Q 
such that 

p= sup \u(x) -u*\ < S , 5 = Tot.Var.{u} <8 , (2.8) 

x€[a,b] 

there exists an e- approximate front tracking solution u = u(t,x) of (1.1), with u(0,x) = u(x), such 
that 

sup |«(3r, x)-u*\< C5 2 , Tot.Var.{u(3r)} < C5 2 . (2.9) 

xE[a,b] 



Proof. On the domain (t,x) G [0,r] x [a, b], we construct u as an e-approximate front tracking 
solution in such a way that, whenever a front hits one of the boundaries x = a or x = b, no reflected 
front is ever created (fig. 3). Since all fronts emerging from the initial data u at time t = exit 
from [a, b] within time r, it is clear that u(t) can contain only fronts of second or higher generation 
order. In other words, the only fronts that can be present in u(t, •) are the new ones, generated 
by interactions at times t > (the dotted lines in fig. 3). Therefore, using the interaction estimate 
(7.69) in [4] we obtain 

sup \u(t, x)-u*\ =0(l)-(p + 8) Tot.Var. {u{t) } = 0(1) ■ 5 2 . (2.10) 

x€ [a,b] 
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figure 3 



We now apply a similar procedure as in the proof of Lemma 1, and construct a solution on the 
interval [r, 3r] in such a way that u(3t) ps u* . More precisely, to construct u on the domain 
[r, 2t] x [a, 6], consider the state u" implicitly defined by (2.2), with t> = u(t, b— ), v' = «*. On 
a forward neighborhood of the point (r, b) we let -U coincide with (a front-tracking approximation 
of) the solution to the Riemann problem 



This procedure will introduce at the point (r,b) a family of wave- fronts of families i = l,...,p, 
whose total strength is 0(1) ■ (p + S). Because of (2.4), all these fronts will exit from the boundary 
x = a within time It. Of course, they can interact with the other fronts present in u(r, ■). In any 
case, the total strength of fronts in n(2r, •) is still estimated as 




if x < b, 
if x > b. 



Tot.Var.{u(2r)} = 0(1) • 5 2 . 



(2.11) 



Next, to define u for t 6 [2r, 3r], consider the state v'" implicitly defined by 



n(2r, o+) = *„(<7„) o • • • o % +1 (<r p+1 )(v'"), 
u* = y p (a p )o...o^ 1 (a 1 )(v'"). 



(2.12) 
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On a forward neighborhood of the point (2r, a) we let u coincide with (a front-tracking approxi- 
mation of) the solution to the Riemann problem 



u(2r,x) = {ffiT,a+) il r > a. 



if x < a. 

This procedure introduces at the point (2r, a) a family of wave-fronts of families i = p + 1, . . . , n, 
whose total strength is 0(1) • (p + 5). Because of (2.4), all these fronts will exit from the boundary 
x = b within time 3r. Of course, they can interact with the other fronts present in n(2r, •). In any 
case, the total strength of fronts in u(3r, •) is still estimated as 

Tot.Var.{u(3r)} = 0(1) • 5 2 . (2.13) 

Moreover, the difference between the values u(3t, x) and u* will be of the same order of the total 
strength of waves in u(r, •), so that the first inequality in (2.9) will also hold. □ 

Proof of Theorem 1. Using the same arguments as in the proof of Lemma 1.1, for every e > 
we can construct an e-approximate front tracking solution u = u(t, x) on [0, 27VY] x [a, b] such 
that 

sup \u(2Nt,x) -u*\ = 0(1) -5, Tot.Var.{u(2iVr)} = 0(1) -5. (2.14) 

x£[a,b] 

Choosing 5 > sufficiently small, we can assume that, in (2.14), 0(1) • 5 < 5$ < 1/C, the constant 
in Lemma 2. Calling T = 2Nt, we can now repeat the construction described in Lemma 2 on each 
interval [T + 3fcr, T + 3(k + l)r] . This yields 

sup \u(T + 3kT,x)-u*\<6 k , Tot.Var.{u(T + 3fcr)} < S k , (2.15) 

x€[a,b] 

where the constants Sk satisfy the inductive relations 

4+i < Cb\. (2.16) 

Choosing a sequence of e-approximate front tracking solutions u e satisfying (2.15)-(2.16) and tak- 
ing the limit as e — > 0, we obtain an entropy weak solution u which still satisfies the same estimates. 
The bounds (1.7)-(1.8) are now a consequence of (2.15)-(2.16), with a suitable choice of the con- 
stants Co, K. □ 



3 - Decay of positive waves 

Throughout the following, we consider a 2 x 2 system of conservation laws 

u t + f(u) x = 0, (3.1) 

satisfying the assumptions (H). Following [6], p. 128, we construct a set of Riemann coordinates 
(wi,W2) ■ One can then choose the right eigenvectors of Df(u) so that 

r,(n) = ^, —± =D \ i . ri>0 i = 1,2. (3.2) 
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It will be convenient to perform most of the analysis on a special class of solutions: piecewise 
Lipschitz functions with finitely many shocks and no compression waves. Due to the geometric 
structure of the system, this set of functions turns out to be positively invariant for the flow 
generated by the hyperbolic system. We first derive several a priori estimate concerning these 
solutions, in particular on the strength and location of the shocks. We then observe that any BV 
solution can be obtained as limit of a sequence of piecewise Lipschitz solutions in our special class. 
Our estimates can thus be extended to general BV solutions. 

Definition 1. We call U the set of all piecewise Lipschitz functions u : M i— > M 2 with finitely 
many jumps, such that: 

(i) at every jump, the corresponding Riemann problem is solved only in terms of shocks (no 
centered rarefactions); 

(ii) no compression waves are present, i.e.: Wi tX (x) > at almost every x G M, i = 1, 2. 

The next lemma establishes the forward invariance of the set U. 

Lemma 3. Consider the 2x2 system of conservation laws (3.1), satisfying the assumptions 
(H). Let u = u(t,x) be the solution to a Cauchy problem, with small total variation, satisfying 
u(0, •) G U. Then 

u(t, ■) G U for all t > 0. (3.3) 

Proof. We have to show that, as time progresses, the total number of shocks does not increase 
and no compression wave is ever formed. This will be the case provided that 

(i) The interaction of two shocks of the same family produces an outgoing shock of the other 
family. 

(ii) The interaction of a shock with an infinitesimal rarefaction wave of the same family produces 
a rarefaction wave in the other family. 

Both of the above conditions can be easily checked by analysing the relative positions of shocks 
and rarefaction curves. We will do this for the first family, leaving the verification of the other case 
to the reader. 

Call a I— ► R\{a) the rarefaction curve through a state uq , parametrized so that 

\i(Ri(a)) =X 1 (u ) + a. 

It is well known that the shock curve through uq has a second order tangency with this rarefaction 
curve. Hence there exists a smooth function c\ (a) such that the point 

St(a) = R^a) + ci(ff)yr2(«o) 

lies on this shock curve, for all a in a neighborhood of zero. From the Rankine-Hugoniot equations 
it now follows 

X (a) = (^f{Ri(o-)+c 1 (a)(a 3 /6)r 2 (u )) - /(«o)) A ^(a) + Cl (a)(a 3 /6) r 2 (u ) - u ) = . (3.4) 
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Differentiating the wedge product (3.4) four times at a = and denoting derivatives with upper 
dots, we obtain 

^(0) = 4[Ai(tt ) Ri (0) + 2/21(0) + A 2 (n ) ci(0)r 2 (uo)] A i?i(0) 

+ 6[Ai(tt )/ii(0) + J Ri(0)] A i?i(0) + 4Ai(u )i?i(0) A [ R\ (0) + c(0)r 2 (n )] 
= 4(A 2 (u ) - Ai(«o))ci(0) r 2 (u ) A n(uo) + 2(Dn • ri)(u ) A n(uo) 
= 0. 

Hence 

Cl (0)= <0 , (3.5) 

u ; 2(A 2 -Ai)(n Ar 2 ) V ' 

By (3.5), the relative position of 1-shock and 1-rarefaction curves is as depicted in fig. 1. By the 
geometry of wave curves, the properties (i) and (ii) are now clear. Figure 4a illustrates the inter- 
action of two 1-shocks, while fig. 4b shows the interaction between a 1-shock and a 1-rarefaction. 
By ui,u m ,u r we denote the left, middle and right states before the interaction, while u' m is the 
middle state after the interaction. In the two cases, the solution of the Riemann problem contains 
a 2-shock and a 2-rarefaction, respectively. □ 




The next lemma shows the decay of positive waves for solutions with small total variation, 
taking values inside U. 

Lemma 4. Let u = u(t,x) be a solution of the Cauchy problem for the 2x2 system (3.1) 
satisfying (H). Assume that 

u(t,-)eU t>0. (3.6) 

Then there exist n, 5 > such that if Tot.Var.(u(t, •)) < 5 for all t, then its Riemann coordinates 
(wi,W2) satisfy 

0<w itX (t,x) <j, t > 0, i = 1,2. (3.7) 

Proof. We consider the case i = 1. Fix any point (t, x). Since centered rarefaction waves 
are not present, there exists a unique 1-characteristic through this point, which we denote as 
1 1-> xi(t; t, x). It is the solution of the Cauchy problem 

±(t) = Ai (u(t, x(t))) , x(t) = x. (3.8) 
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The evolution of W\ tX along this characteristic is described by 

d , N . . 9Ai 2 d\\ 

— Wi >x {t,Xi(t)) = Wl, xt + XlWi >xx = -{XlWl, x ) x + XlW ljXX = -q^ W 1,x - Q^m,xW2,x- 

Since the system is genuinely nonlinear there exists k\ > such that dXi/dwi > k\ > 0, hence 

^w hx {t,xi(t)) < -k x w\ )X + 0(1) • w ltX w 2 , x - (3.9) 

Moreover, at each time t a where the characteristic crosses a 2-shock of strength |<j Q | we have the 
estimate 

w ltX (t a +) < (1 + 0(1) • \<r a \)w ltX (t a -). (3.10) 

Let Q(t) be the total interaction potential at time t (see for example [4], p. 202) and let V^t) be 
the total amount of 2-waves approaching our 1-wave located at x\(t). Repeating the arguments in 
[4], p. 139, we can find a constant Co > such that the quantity 

T(t) = Vi(t) + CoQ(t), t>0, 
is non-increasing. Moreover, for a.e. t one has 

< -|A 2 - Ai||u;2,x|(i,«i(i)) > 
while at times t a where a?i crosses a 2-shock of strength |cr a | there holds 

T(t a -) < T(t a +) - |<7 a |. 

Call = u>i ja; (i,xi(t)). By the previous estimates, from (3.9) and (3.10) it follows form 

W(t) < -hW 2 ^) - CT(t)W(t), (3.9) 

W(t a +) - W(t a -) < C[T(t a +) - T(t a -)]W(t a -), (3.11) 
for a suitable constant C. We now observe that 

p -CT(i) 

y(t) = 



k,e- cr ^ds 



is a distributional solution of the equation 

y = -k iy 2 -Ct(t)y, 

with y(t) — > oo as t — > 0+. A comparison argument now yields W(i) < Since T is positive 
and decreasing, we have 

1 1 e CT (°) 
W(t) < W(t) < , < 



Jo 



for all t > 0. This establishes (3.7) for i = 1, with k = e CT (°) /A^. The case i = 2 is identical. □ 
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We conclude this section by proving a decay estimate for positive waves, valid for general BV 
solutions of the system (3.1). For this purpose, we need to recall some definitions introduced in 
[5]. See also p. 201 in [4]. 

Let u:Mi—> M 2 have bounded variation. By possibly changing the values of u at countably 
many points, we can assume that u is right continuous. The distributional derivative fj, = D x u is 
a vector measure, which can be decomposed into a continuous and an atomic part: fi = (j, c + fj, a . 
For i = 1,2, the scalar measures fi l = [i\ + fi l a are defined as follows. The continuous part of // is 
the Radon measure /i* such that 

J 0d/4 = J <j)h(u) ■ dfi c (3.12) 

for every scalar continuous function cf> with compact support. The atomic part of /j 1 is the measure 
fi l a concentrated on the countable set {x a ; a = 1,2,...} where u has a jump, such that 

a4({ x 4) = <r a ,i = Ei(u(x a -), u(x a +)) (3.13) 

is the size of the i-th. wave in the solution of the corresponding Riemann problem with data u(x a ±). 
We regard fi l as the measure of i-waves in the solution u. It can be decomposed in a positive and 
a negative part, so that 

= tf+ - /J- , |//| =//+ + //". (3.14) 

The decay estimate in (3.7) can now be extended to general BV solutions. Indeed, we show that 
the density of positive i-waves decays as n/t. By meas(J) we denote here the Lebesgue measure 
of a set J. 

Lemma 5. Let u = u(t,x) be a solution of the Cauchy problem for the 2x2 system (3.1) satisfying 
(H). Then there exist n,5 > such that if Tot.Var.(u(t, •)) < 5 for all t, then the measures , 
of positive waves in u(t, •) satisfy 

lA + {J) <^meas(J) (3.15) 
for every Borel set J C M and every t > 0, i = 1, 2. 

Proof. For every BV solution u of (3.1) we can construct a sequence of solutions u u with u v — > u 
as v — > oo and such that u v (t, •) € U for all t. Calling (w", to^) the Riemann coordinates of u„, by 
Lemma 4 we have 

0<wl x (t,x)<j, i > 0, i = 1,2, v>\. (3.16) 

For a fixed t > 0, observe that the map x i— ► (t, x) has upward jumps precisely at the points x a 
where u(t, •) has a 2-shock. Define jl u as the positive, purely atomic measure, concentrated on the 
finitely many points x a where u(t, •) has a 2-shock, such that 

fj-u{{x a }) = Wi(t, x a +) - w\{t, x a -) < C \a a \ 3 (3.17) 

for some constant C. By possibly taking a subsequence, we can assume the existence of a weak limit 
fi v — ^ ft. Because of the estimate in (3.17), the measure ft is purely atomic, and is concentrated 
on the set of points x@ which are limits as v — > oo of a sequence of points x v a where u u (t, •) has 
a 2-shock of uniformly positive strength \a u \ > 5 > 0. Therefore, fi is concentrated on the set of 
points where the limit solution u(t, •) has a 2-shock, and makes no contribution to the positive part 
of n\ + . We thus conclude that the positive part of fi] + is absolutely continuous w.r.t. Lebesgue 
measure, with density < n/t. An analogous argument holds for ji 2+ . □ 
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Corollary 1. Let u = u(t,x) be a solution of the 2x2 system (1-1)- Let the assumptions (H) 
hold. Fixe > and consider the subinterval [a', b'] = [a + e, b — e]. Assume that, at time t = 0, the 
measures /i 2+ of positive waves in u(0, ■) on [a,b] vanish identically. Then, for every t > 
one has 

l4 + (J) <-^-meas(J) (3.18) 
for every Borel set J C [a', b'] and every t > 0, i = 1, 2. 

Indeed, recalling (1.9), the values of u(t, •) restricted to the interval [a' , b'} can be obtained by 
solving a Cauchy problem, with initial data assigned on the whole interval [a, b] at time t — e/X*. 



4 - Proof of Theorem 2 

Lemma 6. In the same setting as Lemma 4, assume that there exists n' > such that 

< w iiX (t,x) < k' te[0,T], i = 1,2. (4.1) 

Let t i ^ x(i) 6e i/te location of a shock, with strength \o~(t)\. There exists a constant < c < 1 
such that 

\a(t)\ > c\a(s)\ , 0<s<t<T. (4.2) 

Proof. To fix the ideas, let u(t, •) have a 1-shock located at x(t), with strength |<r(i)|. Outside 
points of interaction with other shocks, the strength satisfies an inequality of the form 

j t \o-(t) \>-C- (w l!X (t, x(t) + ) + w ltX (t, x(t) - ) w 2 , x (t, x(t) + ) + w 2 , x (t, x(t) - ) ) \a(t) \ . (4.3) 

At times where our 1-shock interacts with other 1-shocks, its strength increases. Moreover, at each 
time t a where our 1-shock interacts with a 2-shock, say of strength |er a |, one has 

\<r(t a +)\ > \a(t a -)\ (1-C'\a a \) . (4.4) 

for some constant C . Assuming that the total variation remains small, the total amount of 2- 
shocks which cross any given 1-shock is uniformly small. Hence, (4.3)-(4.4) together imply (4.2). 

□ 

Lemma 7. Let t i— ► u(t, ■) G U be a solution of the Cauchy problem for a genuinely nonlinear 
2x2 system satisfying (1.11). Assume that there exists n' > such that 

w itX {t,x)<n' te[0,T], i = l,2. (4.5) 

Since no centered rarefactions are present, any two i- characteristics, say x(t) < y(t), can uniquely 
be traced backward up to time t = 0. There exists a constant L > such that 

y{t) - x(t) < L (y(s) - x(s)) 0<s<t<T. (4.6) 

Proof. Consider the case i = 2. By definition, the characteristics are solutions of 

±(t) = A 2 (u(t, x(t))) , y(t) = A 2 (u(t, y(t))) . 
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Since the characteristic speed A2 decreases across 2-shocks, we can write 

*v(.t) 



rv(t) 

y{t)-x{t)<C \wt, x (t,0\ + \w 2 , x (t, 0\ d^ + C Y, ( 47 ) 



where <Si[x, y] denotes the set of all 1-shocks located inside the interval [x(t), y(t)] . Introduce the 
function 

if x < x(t), 



(f)(t,x) 



Moreover, define the functional 



i£x(t)<x<y(t), 
1 if x > y(t). 



*(*) = Yl <t>{t,x a (t)) \<r a (t)\ +C Q(t), 

a€Si 

where the summation now refers to all 1-shocks in u(t, •) and Q is the usual interaction potential. 
Observe that the map t i-» is non-increasing. By (4.5) and (4.7) we can now write 

y(t)-x(t)<C> (1 -*(*)) {y(t)-x(t)) 

for some constant C. This implies (4.6) with L = exp {CT + C"$(0)}. □ 

The next result is the key ingredient toward the proof of Theorem 2. It provides the density 
of the set of interaction points where new shocks are generated. 

Lemma 8. Fix e > and define a" = a + 2e, b" = b — 2e. Consider a 2 x 2 system of 
the form (1.1), satisfying (H). Let u be an entropy weak solution defined on [0,r] x [a, b], with 
r = e/4A*. Let (3.18) hold for all t G [0,r], and assume thatu(0,-) has a dense set of 1-shocks on 
the interval [a" , b"\. Then, for < t < t, the solution u(t, •) has a set of 1-shocks which is dense 
on [a", V — X*t] and a set of 2-shocks which is dense on [a", b"\. 

Proof. By the assumptions of the lemma, there exists a sequence of piecewise Lipschitz solutions 
1 1— ► u v (t) <^U such that u v — > u in L 1 , 

2kX* 

0<< j!e (t,x)<_- i = l,2, u>l, 

and moreover the following holds. For every p > 0, there exists 8 > such that each u u (0, ■) 
(with v large enough) contains at least one 1-shock of strength | ov(0)| > 5 on every subinterval 
J C [a", b"] having length > p. 

To prove the first statement in Lemma 8, fix t € [0, r] and consider any non-trivial interval 
\p,q] C [a", b" — Call s 1— ► s 1— »■ q v (s) the backward characteristics through these points, 

relative to the solution u v . We thus have 

\Pu{s) = \\{U V {S, Pu(s))), { p„(t) =p, 

\ q„{s) = X 1 (u u {s, q v (s))), \ q v {t) = q. 

By Lemma 7, ^(0) — p u (0) > p for some p > independent of v. Hence, each solution u v contains 
a shock of strength |<t„(s)| > 5 located inside the interval [p„(0), g„(0)]. Lemma 5 now yields 
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> cS. By possibly taking a subsequence, we conclude that the limit solution u(t, •) contains 
a 1-shock of positive strength at the point x(t) = lima; I/ (i) G [p, </]. 

To prove the second statement, we will show that the set of points where two 1-shocks in u 
interact and produce a new 2-shock is dense on the triangle 

A = {(t,x); te[0,r], a" < x < b" - X*t}. 

Indeed, let t G [0, r] and p < q be as before. For each v sufficiently large, let t i— > x v (t) be the 
location of a 1-shock in u v , with strength | cr^, (t) | > 5 > 0. Assume x^(-) — »■ x(-) as f — > oo, and 
x v (t) G [p, (?], so that is the location of a 1-shock of the limit solution u, say with strength 
\a(t)\ > 0. 




figure 5 



We claim that the set of times t where some other 1-shock o' impinges on a and generates 
a new 2-shock is dense on [0, i]. To see this, fix < t' < t" < t. For each v sufficiently large, 
consider the backward 1-characteristics y„, z v impinging from the left on the shock x v at times 
t",t' respectively (fig. 5). These provide solutions to the Cauchy problems 



y u (t")=x u (t"), 
z v (t')=x v (t'), 



y v {t) = \i(u v (t,y v {t))), 

z v (t) = Xi(u v (t,z v (t))), 
respectively. Observe that 

*„(0)-y„(0)>p 

for some p > independent of v. Indeed, the genuine nonlinearity of the system implies 



X 1 (u u (t,x u (t)-)) - x v (t) > k u„(t, x„(t)+)) - u„(t, x„(t)-) 



> kS. 



Therefore, 

x v {?) - y„(t') >ff >0, 

for some constant p' > independent of v. By Lemma 6, the interval [y v (0), ^(0)] has uniformly 
positive length. Hence it contains a 1-shock of u u (0, •) with uniformly positive strength |(j„(0)| > 
5 > 0. By Lemma 5, every u v has a 1-shock with strength |<r„(£)| > c5 located along some curve 
1 1— ► x u (t) with 

y v (t) < x u (t) < z„(t) te[0,t'\. 

Clearly, this second 1-shock impinges on the shock x v at some time t v G [t',t"], creating a new 
2-shock with uniformly large strength. Letting v — > oo we obtain the result. □ 
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Proof of Theorem 2. Let Sq > be given. We can then construct an initial condition u(0, •) = 4>, 
with Tot.Var.{0} < So, having a dense set of 1-shocks on the interval [a, b], and no other waves. As 
a consequence, for any e > by Corollary 1 we have the estimate (3.18) on the density of positive 
waves away from the boundary. 

Fix r = e/4A*, and consider again the subinterval [a", b"] = [a + 2e, b — 2e\. We can 
apply Lemma 8 first on the time interval [0,r], obtaining the density of 2-shocks on the region 
[0,t] x [a", b"]. Then, by induction on m, the same argument is repeated on each time interval 
t 6 [mr, (m + l)r] , proving the theorem. □ 

Acknowledgment. The second author warmly thanks professor Benedetto Piccoli for stimulating 
conversations. 
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